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A three-dimensional (3D) highly- directive emission system is proposed to enable beam shaping and beam 
steering capabilities in wideband frequencies. It is composed of an omnidirectional source antenna and 
several 3D gradient-refractive-index (GRIN) lenses. To engineer a broadband impedance match, the design 
method for these 3D lenses is established under the scenario of free-space excitation by using a planar 
printed monopole. For realizations and demonstrations, a kind of GRIN metamaterial is proposed, which is 
constructed by non-uniform fractal geometries. Due to the non-resonant and deep-subwavelength features 
of the fractal elements, the resulting 3D GRIN metamaterial lenses have extra wide bandwidth (3 to 
7.5 GHz), and are capable of manipulating electromagnetic wavefronts accurately, advancing the state of the 
art of available GRIN lenses. The proposal for the versatile highly-directive emissions has been confirmed by 
simulations and measurements, showing that not only the number of beams can be arbitrarily tailored but 
also the beam directions can be steerable. The proposal opens a new way to control broadband 
highly-directive emissions with pre-designed directions, promising great potentials in modern wireless 
communication systems. 



Metamaterials, which contain a mass of artificial subwavelengh unit cells, have sparked exponentially 
increased interest among scientists and engieneers over the past decade due to the fantastic electromag- 
netic (EM) properties that are impossible to be achieved using natural materials. Since the initial 
motivation for scientific pursuit of negative refraction index in 200 1\ physicists and engineers have pioneered 
new ways to guide and manipilate EM waves at will through appropriate choices and arrangements of meta-atoms 
(unit cells), giving rising to the directive^ or enhanced^ emissions, perfect imagings^"^, invisibe cloakings^"^, EM 
tunneling^°, perfect absorbers^ and omnidirectional bending^^, etc. In addition to the exhibited exciting 
phenomena in fundamental science, these conceptually new devices have broken the basic limitations for con- 
ventional devices, and thus provide an initial step towards unique functionalities. Nevertheless, the inherent 
drawbacks of large loss and limited bandwidth due to the resonant profiles of electric or magnetic inclusions as 
well as on occasion complicated fabrication process preclude applications of metamaterials in real world. 

To our relief, gradient-refractive-index (GRIN) metamaterials^^ allow us to arbitrarily tailor the refractive 
index at the unit cell level in the frequencies well below resonances, and thus enable macroscopic EM responses 
with a set of new functionalities^^'^^. In the off-resonance region, the GRIN element exhibits required material 
properties in a wide band with negligible losses, advancing the state of the art of available metamaterials. Recent 
success in GRIN metamaterials has provided extreme controls over EM wavefronts across electrically thin layers, 
and thus facilitates an initial step towards Luneburg lenses, beam collimators, beam bendings, EM wave con- 
centrators, illusion devices, microwave focusing, and highly- directive beam-scanning antennas^^"^^. 

Despite these fruitful achievements, most previous reports were confined to the laboratory demonstration of 
the noted physics, for which the transverse-electric (TE) or transverse-magnetic (TM) polarized EM fields are 
examined extensively in a two-dimensional (2D) parallel-plate waveguide system^^"^l Among them, the coaxial 
probe or point source was commonly employed as an ideal line source, which certainly conflicts with the free- 
space excitation in practice. Moreover, the homogenization of the GRIN materials can be further improved in 
terms of an effective medium. Most importantly, among the rarely reported three-dimensional (3D) lenses^^"^^. 



SCIENTIFIC REPORT; | 4 : 5744 | DOI: 1 0. 1 038/srep05744 



1 




f/GHz CC) f/GHz 



Figure 1 | (a) Topology and (b) retrieved constitutive parameters of the proposed GRIN fractal element. The geometrical parameters of the Sierpinski 
fractal ring are given as di = 7.56 mm, d2 = 0.4 mm, and = 1.16 mm when scale = 1. 



the waveguide was exclusively adopted as the source and thus 
enabled only one desirable directive emission beam. This scheme 
to some extent deviates from the increasing demand on highly pre- 
ferable directive emissions with flexible beam numbers and beam 
directions in a low-cost, low-weight and controllable manner. 
However, this 3D multi-beam highly- directive metamaterial lens 
antenna has not been reported yet so far. 

In this article, the effective excitation and broadband impedance 
match of a 3D GRIN metamaterial lens were comprehensively 
explored under a wave launcher of an actual printed monopole. 
The lens was designed in free space instead of 2D mimicking in 
the parallel-plate waveguide where the electric or magnetic fields 
perpendicular to the 2D plane are constant. A kind of GRIN meta- 
material element constructed based on closed Sierpinski fractal ring 
was proposed, studied, and then employed as a fundamental block to 
build the 3D collimating lens with satisfied index gradients. Due to 
the omnidirectional feature of the monopole, a multi-beam antenna 
system can be readily conceived by employing more lenses. The 
proposed scheme has been verified by a far-field measurement of 
four GRIN lenses embedded with a printed monopole. To the best of 
our knowledge, this is the first realization of a 3D antenna system 
with flexible beam numbers and beam directions for broadband 
highly- directive emissions of more than an octave. 

Results 

Characterization of GRIN fractal metamaterials. The fractal scheme 
employed in this work is from the viewpoint of a much electrically 
smaller GRIN metamaterial element in virtue of the space- filling 
property, which is a new avenue advanced in material design^ 
Similar scheme has been employed to miniaturize the elements of 
composite right/left handed transmission line for a given resonant 
frequency^^'^^. The advantages of using rather a small element can 
be highlighted in twofold. First, the loading of one or two such ele- 
ments may miniaturize considerably the planar circuits or antennas^ 



which is an intriguing characteristic that meets the increasing demand 
on high integration of modern communication systems. Second, the 
volumetric structures made of these periodically arranged elements 
are more feasible to be described by the constitutive parameters under 
the framework of an effective medium^° and the resulting collective 
EM response is less influenced by the parasitic diffraction effects, 
enabling the elements in metamaterials to be analogous to the 
atoms or molecules in natural materials. Therefore, a desirable EM 
functionality can be engineered accurately by tailoring the element 
structures, which would benefit considerably the design accuracy and 
success rate. 

Fig. 1(a) shows the topology of the proposed GRIN metamaterial 
particle and its orientation with respect to the normally incident EM 
waves. The periodicity along x, y and z axis is p^, py, and p^, respect- 
ively. The GRIN element is made of a fractal ring (FR) etched on one 
side of the F4B board (leaving the other side unoccupied) with dielec- 
tric constant 8^ = 2.65, thickness h = 0.3 mm, and loss tangent tana 
= 0.001. Due to the four-fold rotational symmetry of the FR in xoy 
plane, the GRIN element presents a planar isotropy and thus sup- 
ports dual -polarization operation. In this particular design, the 
GRIN element is illuminated by a normally incident plane EM wave 
with polarization in the x direction, and the Sierpinski ring of the 
second-iteration order is utilized by taking a tradeoff between mini- 
aturization and fabrications. Indeed, any zigzag fractal structure with 
strong space-filling property may benefit the compact design. 

In the simulation setup shown in Fig. 1(b), the top-and-bottom 
boundaries of the GRIN element are assigned as perfectly magnetic 
conductors (PMC), whereas the left- and- right boundaries are set as 
perfectly electric conductors (PEC). The impinging electric field 
drives the line inductor in the x direction, and the gap capacitor 
(inter-ring coupling) formed between neighbouring FR structures 
as well as the capacitor formed between close neighbouring sections 
of the FR itself, facilitating an electric resonance modeled by a LC 
resonator which restores and discharges alternately. This can be 
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proved from Fig. SI where the electric resonant frequency 
fQ = l I 2%VlC clearly shifts upwards as and py increase simulta- 
neously and appears to be saturated as p^ and py reach 9 mm. 
Although the decrease of p^ and py ensures increased gap capaci- 
tances, it either affords limited enhancement of capacitance or 
inherently brings about fabrication challenging and tolerances. 
Therefore, it is unadvisable to select this scheme to engineer a com- 
pact element. 

Since the oscillating current forms and flows through the FR struc- 
ture, the abnormal fractal zigzag boundary will succeed in extending 
the current path within a limited volume and thus efficiently lowers 
down the resonant frequency, enabling the element to be less than 
one-tenth of the working wavelength at /q. To afford an efficient 
control over /o in terms of variable inductors and capacitors in a 
simple manner, the Sierpinski ring is progressively scaled up and 
down by manipulating the scaling coefficient within a fixed lattice 
constant. In this regard, the effective parameters below /o change 
accordingly and thus enable us to achieve a controllable arithmetical 
progression of refractive index. 

To illustrate the effects of fractal perturbation, the reflection and 
transmission coefficients of the proposed GRIN fractal particle and 
its conventional counterpart using the square ring (SR)^^ are char- 
acterized by the commercial software CST Microwave Studio. For 
fair comparisons, both rings are with the same scaling coefficient 0.77 
and both elements have the same lattice constant X pyX p^ = 6 X 
6X6 mm^. These S-parameters (not shown for brevity of contents) 
are then utilized to retrieve the effective constitutive parameters 
using the standard retrieval procedure^\ as shown in Fig. 1(c). As 
is observed in both cases, the permittivity curve undergoes an electric 
resonance around/o, whereas the permeability curve suffers an anti- 
resonance. Moreover, an obvious frequency reduction of /o near 
2.9 GHz from 10.25 GHz to 7.34 GHz is inspected when the fractal 
geometry is introduced. Therefore, the GRIN fractal element obtains 
at least 28.3% size reduction for a given operation frequency. 
Moreover, the effective permittivity and permeability in both cases 
exhibit smooth slopes below /o and the permittivity of the fractal 
element suffers a sharper and stronger resonance at/o. As a con- 
sequence, the fractal element features larger index gradient with 
respect to its conventional counterpart. 



Broadband Highly-Directive Emission System. Figs. 2(a)-(d) 
depict the topology and elaborated geometrical parameters of the 
proposed highly- directive emission system or in other words 
integrated lens antenna. As is shown, the antenna is composed of a 
broadband 3D GRIN collimating lens and an omnidirectional 
monopole to launch the lens. The distance from the monopole to 
the lens is D. Distinguished from the curved Fresnel lens which is 
often implemented by polystyrenes, the GRIN lens proposed here 
preserves flat surfaces and contains both core lens and coat lens 
(impedance matching layer) to guarantee a wide operation band. 
Moreover, the GRIN lens is actualized by the proposed fractal 
metamaterials discussed previously, making current implementa- 
tion distinguished from any earlier designs. Here, a planar 
microstrip trapezoid printed monopole deposited on a F4B 
substrate with = 3.5, h = 1.5 mm and tana = 0.001 is explored 
as the source to achieve controllable uniform beams along several 
predicted directions with high integration, low weight, low cost, and 
broad band. To guarantee the bandwidth of the final antenna system, 
both the lens and monopole should operate in a wide band with fully 
overlapped frequencies. The wide bandwidth of the monopole is 
guaranteed by optimizing the geometrical parameters, whereas the 
design principle for a broadband lens is derived in the Methods 
section. Fig. 2(e) depicts the distributions of the theoretically 
calculated refractive index with discrete (solid line) and continuous 
(dashed line) profiles using the derived approach. As can been seen, 
the discrete index of the core lens varies gradually from 2.07 to 1.31, 
while that of the coat lens changes from 1.44 to 1.14. Notice that the 
discretized constitutive parameters have no penalty on lens function 
and performances^^ but enable the lens implementation more 
feasible and practical. 

Although the inhomogeneous index gradient can be easily imple- 
mented by drilling isotropic holes in a dielectric slab^^, the high 
fabrication cost, incompatibility with available printed circuit board 
(PCB) technology, and finite types of available aiguilles with discrete 
aperture dimensions limit the appUcations. Here, the proposed 
GRIN fractal particle is employed to realize the index gradient pre- 
sented in Fig. 2(e) due to the unique EM response discussed in 
previous Section. Since the fractal element meets well the require- 
ment for an effective medium and the material parameters exhibit a 
smooth dispersion in the off- resonant region, the undesirable distor- 
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Figure 2 | Topology and required index gradient of the broadband 3D GRIN lens antenna system, (a) The perspective view of the system, (b) The 
illustration of monopole parameters, (c) The side view and (d) bottom view of the system, (e) The refractive index of the core and coat lens, in which tcore 
= 30 mm, tcoat — 12 mm, S = 90 mm, Li^^s — 160 mm, = 1.5 mm, L = 30 mm, w = 2.6 mm, Hq = 16 mm, hi = 17.5 mm, = 12.5 mm, d= 26 mm, 
and b = 16 mm. 
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Figure 3 | Photographs of the fabricated GRIN slab lens which occupies 
an area of 170 X 170 X 54 mm^. (a) The coat layer and core layer, (b) The 
side view, (c) The perspective view. 

tion of optical quality and significant impedance change between two 
adjacent layers suffered in traditional approach will be mitigated or 
completely eliminated. Moreover, the EM wave outgoing from the 
monopole illuminates the lens with normal incidence along z axis 
and polarization along)/ axis. This EM condition coincides well with 
that presented in Fig. 1(a), and thus will render the desirable electric 
response discussed previously, and in turn the efficient excitation of 
the resulting lens. 

Here, several guidelines are highlighted for an easy design. First, 
the larger the lattice constant along the propagation is, the smaller 
index gradient and lower index values become. Second, the thicker 
the dielectric board is, the larger permittivity while smaller permeab- 
ility are, leading to smaller index and deteriorative impedance 
matching. Finally, the larger is, the larger Ur is. The unit cells of 
the core and coat lens are built on two kinds of F4B boards with (g^, /z, 
tana) = (2.65, 0.5 mm, 0.001) and (2.2, 0.5 mm, 0.001), correspond- 
ing to V12 X V12 X V20 and X^lXl X X^lXl X VlO, respect- 
ively. Here, is the working wavelength in free space at 5 GHz. 

For design and characterization, the GRIN element, the lens, and 
the integrated antenna system are numerically evaluated in CST. The 
retrieved effective material parameters for the core and coat lens at 
5.5 GHz are portrayed in Figs. S2 and S3, respectively. In both cases, 
the dimensions of FR sweep from scale = 0.4 to scale = 0.6 in steps of 
0.001. Therein, we observe that the permittivity and refractive index 
vary continuously, while the permeability almost maintains at unity. 
Therefore, the dimensions of FR reduce gradually from the center to 
the edge of the lens. Moreover, the refractive index of the core lens 
increases progressively from 1.37 to 2.09 while that of the coat lens 
increases from 1.17 to 1.57. These index gradients completely fulfill 
the theoretical requirements for the lens. Most importantly, the ima- 
ginary parts of the permittivity, permeability and refractive index in 
both cases are negligible across the entire sweeping range, indicating 
a low loss. 

Given the parametric sweeping results, the lens can be exactly 
constructed based on the geometrical mapping process using a 
root- finding algorithm. Fig. 3 depicts the photograph of the fabricated 
3D GRIN lens which is composed of ten core layers sandwiched by 
two coat layers on each side. In the radial direction, each concentric 
region shown in Fig. 2(d) is reaUzed by two GRIN elements with 
identical geometrical parameters. Hence a total of 32 X 32 metallic 
patterns are positioned on either side of each F4B dielectric board. An 
additional 2.5-mm-thick foam plate with near-unity permeability and 




(b) 





Figure 4 | The scheme and prototype of the proposed multi-beam highly- 
directive emission systems, (a) The perspective view and (b) the side-view 
of a two-beam system, (c) The perspective view of a four-beam system, (d) 
The side-view of the four-beam systems, in which the four lenses are 
rotated by 0 = 0°, 0 = 30°, Q = 45°, and 9 = 60°, respectively, (e) The 
prototype of the four-beam highly-directive emission system. 



permittivity is placed between every two neighboring core layers to 
guarantee pz = 3 mm, whereas an additional 5.5-mm-thick foam 
plate is inserted between every two coat layers to guarantee = 
6 mm. By tightly stacking these dielectric boards and foam plates 
alternatively and reinforcing them through adhesives, a 3D GRIN slab 
lens occupying an area of 170 X 170 X 54 mm^ is implemented. The 
lens is packaged in hand-made foam to provide supporting frame and 
can be made more compact when smaller index gradients or less 
concentric layers are used. However, this would either degrade the 
accuracy of lens implementation or deteriorate the lens performances. 

The single-beam highly- directive emission system can be directly 
generalized to a multi-beam one. Fig. 4 illustrates the schematic of 
various multi-beam emission systems and the fabricated prototype of 
a four-beam lens antenna for demonstration. As can be seen, two or 
four GRIN lenses are arranged around the source monopole with 
identical D. The highly- directive emissions can be controlled in fol- 
lowing four manners. First, the number of beams can be engineered 
arbitrarily by appropriately distributing more GRIN lenses. Second, 
the beam directions can be arbitrarily manipulated by tailoring the 
included angles between neighboring lenses. Third, the beams can be 
engineered rotated by placing the lenses on a mechanically rotated 
platform isochronous with a time sequence and fixing the feeding 
antenna. Finally, the beam can be tailored nonuniform by varying D 
and index gradient of the lenses. Here, the direction of beam is 
defined corresponding to the orientation angle 0 with respect to z 
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Figure 5 | Simulated 3D far-field radiation patterns at 5.5 GHz for (a) conventional monopole, (b) single-beam, (c) two-beam, and (d-g) four-beam lens 
antennas with 6 = 0° (d), 6 = 30° (e), 6 = 45° (f), and 6 = 60° (g), respectively, (h) Reflection coefficients of these antennas. 



Numerical and experiemtal results. To validate the design method 
for the lens antenna, the gain is plotted as a function of D, see Fig. S4. 
The maximum gain of 14 dB occurs around D = 90 mm which is 
used to design the index gradient. Consequently, the rationality of 
approximating the monopole as a point source and effective design of 
the lens are confirmed. Higher gain of more than 19.3 dB and larger 
aperture efficiency are achieved at 5 GHz when a horn antenna with 
length of D = 90 mm and aperture of the lens is employed as the 
source at the cost of multiple beams. This level of gain approximates 
to the maximum gain of 20 dB according to Gmax = 47iL^^ /y^o for 
which D should be sufficiently long. Therefore, the relative low gain 
of the lens antenna is attributable to the utilized omnidirectional 
monopole. In comparison with multi-beam emission system using 
several aperture antennas or arrays, our approach features moderate 



antenna gain, low loss, broad band, and requires only a single feed 
without using complicated feeding networks such as power dividers. 
Moreover, the distance D in our scheme can be significantly reduced 
relative to the bare horn for the same level of uniform phase 
distribution, enabling a compact multi-beam antenna. 

To demonstrate the phase-coUimating capability of the GRIN lens, 
the aperture fields in xoy plane for antennas with and without the GRIN 
lens are calculated and the electric field Ey component is plotted, see Fig. 
S5. Clearly, desirable phase-correction fimctionality can be identified 
from the uniform phase fi-onts across the entire aperture except for the 
edge of the lens. On the contrary, the phase fi-ont suffers significant 
fiuctuations in the bare monopole case. Notice that the variation of 
phase front on the lens edge poses little infiuence on the radiation 
performances since most energy radiates through the central part. 
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Figure 6 | Simulated far-field radiation patterns and gain of the proposed highly-directive emission system under different scenarios. The snapshots 
are plotted at different frequencies, (a) Single beam, (b) Two beams. Four beams with (c) 6 = 0°, (d) 6 = 30°, (e) 6 = 45° and {{) 6 = 60°. (g) The 
frequency- dependent gain of these antennas. 



Fig. 5 illustrates the simulated 3D far-field radiation patterns 
and reflection coefficients Sn of the monopole with and without 
GRIN lenses, aiming to provide an intuitive view on the controllable 
highly- directive emissions. Different scenarios of single beam, 
two beams, and four beams are extensively examined with different 
directions at 5.5 GHz. As is expected from Figs. 5(a)-(g), the 
omnidirectional patterns of the monopole in xoz plane (H plane) 
are transformed to several narrow-beam patterns with high direc- 
tivity. With respect to the bare printed monopole, an average of 
10 dB directivity enhancement is achieved along lens directions 
in all cases, which is a very attractive and challenging feature. 
Moreover, the monopolar pattern with null radiation at broadside 
still holds, indicating that these lenses only rearrange the energy 
distribution in H plane. Most importantly, it is demonstrated for 
the first time that the number of beams in free space can be tailored 
at will and the beams can be arbitrarily directed by simply steering 
the orientations of the lenses from an embedded actual printed 
monopole. The multi-beam directive emissions deviate from the 
multi-beam generation in^^'^^'^^ which are still constrained to sur- 
round a ideal probe with zero-refractive-index or GRIN metamater- 
ials in a 2D waveguide. Our multi-beam control scheme provides 
more versatility, flexibility and stronger signal intensity, while also 
maintains easy implementation in practice. From Fig. 5(h), it is 
learned that the lenses introduced in all cases do not deteriorate 
the broadband impedance-matching performance of the printed 
monopole, for which the Sn is better than —10 dB from 3.5 to 
7.5 GHz. 



Fig. 6 gives the simulated far- field patterns in H plane and the 
realized gain at a number of frequencies from 3 to 8 GHz. Here, the 
antenna gain in all cases is calculated by taking losses into account. 
From Figs. 6(a)-(f), we learn that the sidelobes increase appreciably 
at the lower and upper edges of the passband, which is attributable to 
the dispersive GRIN material parameters. In this scenario, the con- 
stitutive parameters of the lens deviate slightly from the theoretically 
calculated ones when the frequencies deviate from the centre fre- 
quency 5.5 GHz at which the lens is exactly designed. Nevertheless, 
the slightly deteriorative sidelobe in the off-center band does not 
preclude the system from real-world applications since the directivity 
of main beams does not decay until 6.5 GHz. From Fig. 6 (g), an 
obvious gain enhancement is observed across the whole band in all 
cases and the minimum value is found to be 6.5 dB at 3 GHz. The 
gain decay at the lower frequency end is due to the size of the aperture 
in terms of the wavelength decreases significantly at that end, 
whereas the dispersive GRIN material parameters account for the 
raised return loss and material losses and thus give rise to the gain 
decay at the upper frequency end. Despite this, the radiation effi- 
ciency of the lens antenna r\ calculated as the proportion of the gain 
and directivity is larger than 90%, further demonstrating that the 
losses induced by the lens are extremely small. Moreover, the gain 
is very appealing and comparable with respect to aforementioned 
bare horn antenna, see the dashed blue curve. This is especially true 
for the gain at high frequencies. 

To provide physical insight regarding the highly- directive emis- 
sions, the near-field mapping is performed in 3D free-space mea- 
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Figure 7 | Measured electric-field distributions and reflection coefficients. The snapshots are plotted in free space behind the bare monopole and lens in 
steps of 2 mm at several selected frequencies. 
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Figure 8 | Simulated and measured radiation patterns for the (a) -(c) single-beam and (d)-(f) four-beam antenna at 3.5 GHz (the left column), 
5.5 GHz (the middle column), and 7.5 GHz (the right column). Note that the results of the bare monopole (without the lens) are given for comparisons. 



surement apparatus. In the experimental setup, see Fig. S6, the 
broadband monopole placed on an aluminum plate is 90 mm away 
from the GRIN lens and is embedded in foam to provide mechanical 
support. It functions as a transmitter and moves together with an 
electronic step motor automatically along two dimensions by scan- 
ning an area of 200 X 200 mm^. At the front side of the lens, a 
1.19 mm-diameter semi-rigid coaxial probe is fixed stationary to a 
ferric bracket to detect the electric fields. Both transmitting and 
receiving antennas are connected to the two ports of a phase-sens- 
itive network analyzer (Agilent PNA-LN5230C) via two coaxial 
cables and are placed at the waist of the GRIN lens to avoid trun- 
cation effects of finite lens dimension. 

Fig. 7 plots the reflection coefficients and measured free-space 
electric-field distributions. The successfiil conversion of quasi- spher- 
ical waves to plane waves is clearly observed from 4.5 to 8 GHz in 
steps of 0.5 GHz by comparing the wavefronts outgoing from the lens 
antenna and bare monopole, respectively. Here, the results at the 
lower edge of passband are not given due to the large wavelength at 
these frequencies. In this regard, the long wavefront exceeds the scan- 
ning capability of available near-field measurement apparatus, mak- 
ing it impossible to clearly identify a complete wavefront. 
Nevertheless, the highly-directive emissions at these frequencies can 
be evidenced from the measured far-field patterns. Since the GRIN 
elements operate far from their resonant frequencies, the loss induced 
by the sample is extremely low and is hardly indicated from the 
measurements. Therefore, almost identical intensity can be inspected 
from the color maps of outgoing waves in all cases. The nonideal 
monopole with finite ground dimensions gives rise to the uneven 
and unsmooth wavefronts. Moreover, the measured reflection coeffi- 
cients are in reasonable agreement with simulations, revealing that the 
return loss of the bare monopole and integrated lens antennas under 
different scenarios is better than — 10 dB from 3.2 to 7.75 GHz, which 
is more than an octave. Therefore, the accuracy of employing GRIN 
fractal elements in constructing bulk GRIN lenses is verified. 



To experimentally validate the appealing far-field characteristics of 
the 3D antenna system. Fig. 8 compares the simulated and measured 
H-plane radiation patterns of both single-beam and four-beam anten- 
nas from 3.5 to 7.5 GHz in steps of 2 GHz. The radiation patterns in 
other frequencies can be observed from Fig. S7. All patterns are 
normalized to the co -polarization of the designed monopole. As is 
inspected, the results between simulations and measurements are in 
reasonable consistency. In both single-beam and four-beam cases, the 
gain enhancement is observed more than 10 dB along certain orien- 
tations. The measured cross-polarization radiations are at least 15 dB 
lower than the co -polarization ones at most frequencies and become 
worse at both edges of the passband. The nonuniform intensity 
(amplitude fluctuations) of the four main beams is attributed to the 
slightiy broken uniformity of the monopole pattern which is induced 
by the insufficiently large ground plane. The appreciably higher cross- 
polarization radiations in experiments, and slight discrepancies of co- 
polarization patterns refiection coefficients (especially the shift of 
reflection zeros) between simulations and measurements are partially 
attributable to the tolerances that are inherent in fabrications, and 
partially to effects of supporting foam and adhesives that are not 
considered in simulations. Moreover, the misalignment of transmit- 
ting and receiving antennas and the nonideal background envir- 
onment also account for the observed deviations. Nevertheless, 
these tolerances pose litfle influence on lens function. 

Discussion 

To conclude, we have succeeded in demonstrating a highly- directive 
emission system which possesses the merits of broad band, low loss, 
high gain, and flexible beam numbers and beam directions. 
Numerical and experimental results are in good consistency, show- 
ing that the emission system exhibits a near 10 dB gain enhance- 
ment, and moderate cross-polarized radiations which are 15 dB 
lower than the co-polarized radiations over an octave in super- 
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extended C band. Therefore, the effective design for a broadband 3D 
slab lenses under the free-space excitation of a planar printed mono- 
pole and the advantages of using fractal strategy in 3D GRIN lens 
design without increasing the manufacturing complexity have been 
confirmed. The proposed scheme introduces an alternative avenue to 
improve the gain of microstrip antennas from a material perspective. 
Other omnidirectional source with uniform patterns can be 
employed for further improvements. 

Methods 

Suppose that the aperture of the monopole is negligible with respect to that of the 3D 
lens, the monopole can be approximated as a point source in design and the dif- 
fraction in terms of truncation effects can be eliminated. Despite this, the design 
process is more complicated than 2D design by using ideal point source, which will be 
discussed later. Denote the thickness and side length of the lens as t = fcore + 2tcoat 
and Liens> respectively and assume that the wave inside the lens travels parallel to the 
optical axis without a lean. Here, tcore and tcoat are thicknesses of the core and coat lens, 
respectively. To guarantee the phase-collimating performance of the lens, the optical 
path from the source to any point in the upper surface of the GRIN lens should be the 
same, yielding the same phase delay. 

D ^ not _ Vr^ + S^ ^ rirt 

where Hq is the refractive index at the center of the lens, and Ur is the radially 
distributed index as a function of radius r. Then, rir can be immediately obtained from 
Eq. (1). 

nr = no + (D-Vr^TD^)/t (2) 

Although the required index of the lens is obtained, it is impossible to achieve the 
rigorous criterion 8r = fir= rir in broadband and there is a large index fluctuation at 
the interface of lens and air, resulting in a severe impedance mismatch. To relax the 
fluctuation and thus enable a smooth transition of index and desirable impedance 
matching, the aforementioned one-layer lens is substituted by a core lens sandwiched 
by two identical coat lenses^^ each functioning as a one- stage X/A impedance trans- 
former, see Fig. 2(c). To guarantee simultaneously the beam-collimating function, the 
optical path for the three-layer lens should meet the following criterion. 

"corefcore + 2ttcoatW = ^r^ (3) 

Since the effective permeability in this particular design maintains as unity across the 
entire band, we have ricore = y/^cove and ricoat = A/Scoat- Then Scoat = A/ScoreSo and 
^coat = \/^We^ can be immediately obtained by taking the one- stage X/A transformer 
condition Zcoat = yZ^meZo- By inserting these equations into Eq. (3), the radially 
distributed index of the coat and core lens can be achieved after some simplifications. 

2W (4) 

2 

^core ^coat ^coat 

Here, the thickness of the coat lens can be roughly determined corresponding to Xq/A, 
whereas that of the core lens should be designed as Xq/I to guarantee desirable phase- 
collimating performances. In actual realizations, the 3D collimating lens is discretized 
into eight concentric circular regions each having a specific refractive index. Hq is 
cautiously chosen as Hq = 1.8 to ease the fabrications in achieving large index gra- 
dients. Since the monopole is not an ideal point source, it is hard to strictly align the 
radiation center of the monopole with that of the lens, making the simple ray tracing 
do not always hold. To make the approximation more accurately, the aperture of the 
lens and D should be sufficiently large and should fulfill the condition 

^Jd'^+LI^Ja - + {Liens - Lf /"^ < (7 and ^D^-\-L^/4: -D<g. The parameter 
o defines the accuracy criterion and is chosen as (7 = XqI6 in this particular design. 
Moreover, the relative position (xoy plane) of the monopole relative to the lens is 
optimized to have a maximum peak gain for a given D. Here, D plays a crucial role in 
determining the antenna performances such as gain when the position of the 
monopole is determined. To afford efficient excitation and maximum aperture effi- 
ciency, it is designed and optimized as D = 90 mm by comparing the performances of 
antennas with different D and index gradients. In this connection, the radiation center 
of the outgoing wave can be tailored at the geometrical center of the lens, which 
guarantees a desirable gain and directivity. 
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